Strong-disorder renormalization for interacting non-Abelian anyon systems in two dimensions 
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We consider the effect of quenched spatial disorder on systems of interacting, pinned non-Abelian anyons as 
might arise in disordered Hall samples at filling fractions v — 5/2 or v — Y2l§. In one spatial dimension, such 
disordered anyon models have previously been shown to exhibit a hierarchy of infinite randomness phases. Here, 
we address systems in two spatial dimensions and report on the behavior of Ising and Fibonacci anyons under the 
numerical strong-disorder renormalization group (SDRG). In order to manage the topology-dependent interac- 
tions generated during the flow, we introduce a planar approximation to the SDRG treatment. We characterize 
this planar approximation by studying the flow of disordered hard-core bosons and the transverse field Ising 
model, where it successfully reproduces the known infinite randomness critical point with exponent tp w 0.43. 
Our main conclusion for disordered anyon models in two spatial dimensions is that systems of Ising anyons as 
well as systems of Fibonacci anyons do not realize infinite randomness phases, but flow back to weaker disorder 
under the numerical SDRG treatment. 

PACS numbers: 73.43.-f, 75.10.Nr, 72.15.Rn, 03.65.Vf 



I. INTRODUCTION 



Since von Klitzing's seminal discovery of the quantized 
Hall effect 1 , non- symmetry breaking topological order has 
become an essential part of our understanding of low- 
temperature electronic sy stems P The fundamental feature of 
many of these topological phases of matter is the presence 
of anyonic quasiparticles, whose adiabatic exchange entails a 
nontrivial operation on the state of the system in contrast to 
the signs accumulated by conventional fermions and bosons. 
In the most general case, the exchange of such anyons in- 
duces non-Abelian transformations on a (topologically) de- 
generate manifold of states of the system. There are sev- 
eral candidate systems currently under intense experimental 
scrutiny, which on theoretical grounds have been proposed to 
exhibit the simplest incarnation of such non-Abelian quasipar- 
ticles, so-called Ising anyons. These include the Moore-Read 
state 3 proposed for the fractional quantum Hall liquid at fill- 
ing fraction v — 5/2, p^+ip^ superconductors p heterostruc- 
tures involving topological band insulators 5 6 and certain frus- 
trated magnets.^ An incarnation of slightly more complicated 
anyons, so-called Fibonacci anyons, have been proposed in a 
theoretical description of the FQH state at filling v = 12/5 
based on the Read-Rezayi state. 8 In real samples, however, 
unavoidable impurities pin these particles randomly in space 
and residual microscopic interactions split the various fusion 
channels for the associated topological charge. The collec- 
tive behavior due to this disorder pinning may strongly influ- 
ence thermodynamic transport and topological interferometry 
experiments. It may also impede the use of such systems as 
quantum computers in proposed schemes of topological quan- 
tum computation^. 

The pinned anyon problem in a disordered Hall bar is remi- 



niscent of a quantum spin glass 10 . However, within the Strong 
Disorder Renormalization Group (SDRG) method 1 -^ that we 
use in this article, the non-Abelian character of the fusion 
rules imply that the interactions generated by the SDRG are 
intrinsically topology-dependent. The action of this renor- 
malization group (RG) on fusion in one piece of the sample 
may influence the fusion of particles, and thus their interac- 
tions, elsewhere in the sample, even if they are a priori not 
connected by direct energetic interactions. Hence, many of 
the properties of the SDRG known from applications to tradi- 
tional spin systems, which have simple tensor product Hilbert 
spaces, must be revisited for non-Abelian systems. The one- 
dimensional incarnations of these interacting anyon problems 
have been solved both in the pure 12 and, within the SDRG, in 
the disordered^ES case. In the latter case, the SDRG anal- 
ysis found flows to a hierarchy of infinite randomness fixed 
points whose specific character depends on the underlying 
anyon theory. In this report, we use the SDRG to numeri- 
cally study disordered non-Abelian anyon systems in two spa- 
tial dimensions. We use several different approximations to 
handle the explicitly topology-dependent interactions gener- 
ated by the SDRG renormalization scheme. Our main con- 
clusion from such an SDRG treatment is that in the presence 
of disorder two-dimensional systems of interacting Ising or 
Fibonacci anyons do not realize infinite randomness phases, 
but that these systems flow back to weaker disorder under the 

sdrgP 

In the absence of disorder, these 2-D interacting anyon 
models have been studied previously, and argued to ex- 
hibit gapped topological phases on their owrP^CU. Cou- 
pled with the results of the present study, one might suspect 
that quenched disorder could simply be irrelevant (in the RG 
sense) for two-dimensional anyon models, with all RG flows 
returning to the clean fixed points. Two observations however 



suggest that such a scenario cannot represent the complete pic- 
ture. First, in one dimension, both pure and disordered anyon 
chains exhibit gapless critical phases which are distinct from 
the clean systems and are protected by the topological nature 
of the Hilbert space. Second, by recasting the disordered two- 
dimensional Ising anyon problem in terms of non-interacting 
Majorana fermion zero-modes, some of us^ have recently 
established the presence of a disorder-induced (but not infi- 
nite randomness) thermal metal phase in the phase diagram 
of the pinned disordered 2-D Ising anyon problem. Whether 
an analog of this thermally conducting 2-D phase might also 
be found for the more complicated Fibonacci anyons (e.g. in 
the v — Yljh quantum Hall state) - or any higher level anyon 
model - is an intriguing open question. 

In this report, we first briefly review the physical picture of 
anyons pinned within a two-dimensional quantum Hall sam- 
ple and motivate the effective Hamiltonian describing their 
interactions in Sec. II. In the following section III, we dis- 
cuss the ingredients of a fully two-dimensional strong disor- 
der renormalization group analysis and introduce the planar 
approximation, which we invoke to study the infinite random- 
ness behavior of the Fibonacci anyons. We characterize the 
planar SDRG by comparing its behavior on the well studied 
random transverse field Ising model (TFIM) and hard-core bo- 
son hopping problem and uncover a pathology which we ar- 
gue reflects the flow of the more exotic topological models to 
weaker disorder. We summarize these arguments at the end of 
Sec. Ill Section [TV] provides a more detailed introduction to 
each of the models, the strong-disorder rules and a brief sum- 
mary of the numerical results. Section [V] discusses a few of 
the implementation details for our numerical study. 



II. DISORDERED PINNED ANYON MODELS 

In this section, we provide a heuristic introduction to the 
physics underlying disordered pinned anyon models using the 
v = 5/2 state as an example. For more details, see Refs. 12, 
HBandEl 

Let us consider a Hall bar in a strong magnetic field at 
filling factor v = 5/2 and assume that the fractional quan- 
tum Hall phase associated with this filling is indeed described 
by the Moore-Read Pfaffian state. 3 In an ideal sample tuned 
precisely to v = 5/2, the system forms a uniform electron 
fluid with a gap to quasiparticle (qp) excitations with elec- 
tromagnetic charge e/4 and non-Abelian braiding statistics. 
These quasiparticles are also called Ising anyons or "sigmas", 
in varying contexts. 

Detuning the filling fraction away from the center of the 
5/2 plateau or the introduction of a random electrostatic back- 
ground (eg. due to sample impurities) introduces a finite den- 
sity of qps into the ground state of the system. In the pure, de- 
tuned system, the dilute gas of charged quasiparticles Wigner 
crystalizes into a triangular lattice; in a weakly disordered 
sample, the lattice sites randomly shift toward wells of the po- 
tential. In either case, the orbital (charge) degree of freedom 
of these particles ought to gap out of the low energy spectrum, 
see Figure [T] 



Due to the non-Abelian statistics of the quasiparticles, how- 
ever, this is not the end of the low energy description of the 
system. Assuming the N pinned qps are sufficiently far apart, 
then there remains a manifold of (nearly) degenerate ground 
states that grows exponentially with the number of qps N. 
This degeneracy is in many ways analogous to that of a sys- 
tem of noninteracting spin- 1/2 quasiparticles, which would 
provide the ground state a2 iV spin degeneracy. For the non- 
Abelian qps, the degeneracy depends on the underlying any- 
onic theory and the so-called quantum dimension of the non- 
Abelian degree of freedom in this theory: for Ising anyons it 

asymptotically grows as y 2 , while for Fibonacci anyons it 
grows as (j) N where <j) is the golden ratio <j) = (y/h + l)/2 
(hence the name Fibonacci anyon). These degeneracies are 
not associated with any local observables, but rather with the 
braiding history of the qps. 

We can understand the construction of the topological 
Hilbert space from its fusion rules. For the v = 5/2 phase, 
there are three topological charges: the vacuum 1, the elemen- 
tary quasiparticle a and the fermion ip. That is, measuring the 
net topological charge of any collection of quasiparticles will 
produce one of these three results. The topological charge of 
a single qp is a. 

Thus, we can build up a basis for the topological Hilbert 
space by considering the state space built by successively fus- 
ing together each of the N a 's in the system. The fusion rule 
for a pair of a particles is 



a x a = 1 + if; . 



(1) 



while both the 1 and ip particles act like the identity when 
they fuse with a a. An orthogonal basis for the Hilbert space 
is therefore given by the labelings of a fusion tree (see Fig. [2]) 
that are consistent with the above fusion rules. Notice that this 
description of the Hilbert space requires choosing an ordering 
of the quasiparticles which is implicit in the depiction in the 
form of a fusion tree in Fig. [2] 

- other fusion orderings provide alternative bases which 
may be related by unitary transformations built out of so- 
called 'F-moves' illustrated in Fig. [2] and braid moves (ex- 
changes). 

The above description of the topological Hilbert space cor- 
responds to a collection of qps that are arbitrarily far apart 
- at finite separations, the degenerate manifold of states will 
split in some non-trivial way. This effect can already be seen 
for a pair of anyons where topological charge tunneling 22 will 
result in a splitting of the two possible fusion outcomes ([T]) 
for a pair of sigmas. Generalizing this pair- wise splitting into 
an effective many-body Hamiltonian we arrive at an anyonic 
version 12 of the conventional Heisenberg model given by 



(ij) 



(2) 



where lif- projects onto the ip fusion channel of the ij' th 
pair of cr's and J^ may take any sign and strength depending 
on which channel is locally favored and by how much. The 
tunneling-mediated exchange couplings J^ fall off exponen- 
tially with the interparticle separation l^ relative to the length 




FIG. 1 . (left) Schematic depiction of quasiparticles of size £ localized on a randomly displaced triangular lattice within a two-dimensional 
sample, (right) Schematic spectrum with a bulk gap A and a topological degeneracy split at an exponentially smaller scale in the interparticle 
separation /. 





FIG. 2. (left) Fusion tree for 5 a particles with all possible label- 
ings of each intermediate fusion channel, consistent with the fusion 
of the charges from top to bottom. The four available labelings pro- 
vide a basis for the Hilbert space, (right) Alternative fusion ordering 
providing a different basis for the same space. The bases may be re- 
lated by a unitary transformation given by F-moves and braid moves 
(exchanges). 



scale of the qp wavefunction - roughly the magnetic length. 23 
Thus, all of the physics described herein ultimately lies in the 
narrow band below the gap indicated in Fig.[T] 

It is difficult to get a solid handle on the microscopic 
physics governing J (I). Nonetheless, calculations within sev- 
eral theoretical frameworks have been done p3EI They all find 
an exponentially decaying envelope expected for a tunneling 
mediated process in a gapped quantum liquid, within which 
the favored fusion channel and thus the sign of the J (I) os- 
cillate - akin to an RKKY interaction. Thus, unless the sys- 
tem is an ideal triangular lattice with all kj equal, one expects 
that the Jij indeed have a very broad distribution and strong 
sign disorder. It is reasonable to expect the physics to be cap- 
tured by passing to a disordered ensemble of any on lattice 
Hamiltonians H with independent identically distributed J^ , 
a la the Edwards- Anderson approach to disordered magnets. 
This is also what motivated us to look into strong-disorder ap- 
proaches for a description of the system. 

Any on projector Hamiltonians like ^ depend on geometry 
in a rather subtle way, unfamiliar from spin models. Indeed, 
the notation Ilfj is ambiguous - in general, the pair-wise in- 
teraction of two distinct anyons i and j should be labeled by 
a path connecting the two anyons which indicates on which 
side of the other anyons in the system the interaction is me- 



diated. From the point of view of the fusion basis described 
above, this corresponds to fixing the ambiguity about how to 
change the fusion ordering in order to implement the inter- 
action of two particles not adjacent in the original ordering. 
Within the k — 2 (Ising) anyon theory which we have dis- 
cussed so far, the non-locality does not play an important role 
in our further analysis, but for higher k, it can lead apparently 
uncoupled clusters of anyons to generate explicit couplings 
perturbatively. We will return to this feature in the discussion 
of the Fibonacci rules in Sec. lIVDl 



III. STRONG DISORDER TREATMENT 

We consider the behavior of pinned anyon lattice models 
in the strongly disordered regime. In one spatial dimension, 
such models exhibit a hierarchy of infinite randomness phases 
under the strong-disorder renormalization group (SDRG) for 
SU(2)k anyons indexed by the level /cP^It is natural to ask 
whether such infinite randomness behavior carries over to the 
two dimensional pinned lattice models that motivated their 
study. This would be an especially intriguing result as very 
few two dimensional models are known to flow to such infi- 
nite randomness fixed points. In fact there are only two known 
examples: the transverse field Ising model (TFIM) with ran- 
dom fields and random bonds, which has an infinite random- 
ness critical point (IRCP) separating random ferromagnetic 
and paramagnetic phases 26 ; and the bipartite Majorana hop- 
ping model, which has a marginally stable infinite randomness 
phase. 27 Meanwhile, such simple models as hard-core boson 
hopping flow away from infinite randomness towards weaker 
disorder in two dimensions. 

In one dimension, many infinite randomness fixed points 
(IRFPs) may be found analytically because the strong- 
disorder rules preserve chain geometry and introduce only 
trivial correlations into the coupling distributions P3 This 
holds even for anyonic chains because fusion need not sig- 
nificantly reorder the anyons in the chain. 14 15 In two dimen- 
sions, however, two complications arise: first, the RG rules for 
any model generate lots of next-neighbor bonds which quickly 



render the geometry of an initial 2-D lattice unrecognizable. 
Moreover, these bonds have significant correlations between 
their strength and geometric significance, even as the geome- 
try becomes more obscure. Second, for the anyon models at 
k > 2, the topology of bonds which cross becomes important 
for the generation of renormalized bonds. As the mesh renor- 
malizes, keeping track of such crossings becomes more and 
more problematic. 

The loss of the lattice geometry renders direct analytic treat- 
ments in two dimensions intractable but may be dealt with by 
numerical simulation of the SDRG flow on large instances. In 
this approach, we specify the perturbative rules for integrat- 
ing out strong bonds in a particular model and then iteratively 
apply these rules to decimate large samples while monitor- 
ing the flow of their couplings and geometry. In practice, we 
also need to control the growth of memory requirements by 
some truncation of generated weak couplings. This RG pre- 
dicts its own success or failure: if after an initial transient, we 
discover that the width of the distribution of log couplings 
Pt(\n(Q/J)) falls onto a scaling distribution with a width 
w ~ N~^/ d , then we have found an IRFP with scaling expo- 
nent i/j (here, N is the number of sites remaining in the system 
and d the spatial dimension). If on the other hand, the width 
w shrinks, then the system is trying to flow back to weaker 
disorder and the strong-disorder approach is suspect. 

The effect of crossings on the interaction renormalization, 
which arises for Fibonacci anyons (k = 3), is much more 
problematic for a numerical treatment. Rather than attempt 
to keep track of the topology of crossing bonds, we invoke 
the following approximation: since the original lattice is 2-D 
and we expect the physics to be dominated by nearest neigh- 
bor interactions, we require all renormalized interactions to be 
planar. Thus, on a given renormalization step, newly gener- 
ated bonds are added into the model in order from strongest 
to weakest so long as they do not make the lattice non-planar. 
Clearly, this truncation rule is an approximation to the SDRG 
which is not itself perturbatively controlled. We therefore 
compare the behavior of the SDRG with and without this pla- 
nar approximation on a number of (non-)topological models 
in order to determine if it deems reliable for the case of Fi- 
bonacci anyons. 

Table I summarizes the qualitative behavior of each of the 
models that we have studied using our numerical implementa- 
tion of the SDRG. In order of increasing complexity, these are 
hard-core boson hopping (aka. XX model), Ising anyons in- 
teracting on bipartite and non-bipartite lattices, the transverse 
field Ising model and the pinned Fibonacci anyon model. For 
all of these models except the Fibonacci anyons, we have stud- 
ied the SDRG flow with the planar approximation and without 
in order to understand better the behavior of the approxima- 
tion. The background, RG rules, and quantitative results for 



each of the models appear in more detail in Sec. IV] General 
implementation details will shortly be mentioned in SecfV] 

The primary results of our investigation into the planar ap- 
proximation are as follows: 

• The planar approximation appears to be asymptotically 
valid for the one two-dimensional model that is known 
to flow strongly to infinite randomness — the TFIM at 



its infinite randomness critical point. In particular, the 
existence and scaling properties of the critical point are 
identical with and without the approximation to within 
numerical error. 

• The models which flow back to weaker disorder under 
the usual SDRG are 'stabilized' by the planar trunca- 
tions, in the sense that they find marginal scaling fixed 
points for their bond strength distributions that are inde- 
pendent of system size. In particular, this behavior can 
be seen in the hard-core boson model. 

• These "marginal" fixed points are, however, physically 
spurious. At these fixed points, the distribution of 
strong bond lengths, as measured by the bare geome- 
try of the original lattice positions, saturates the system 
- indeed, the strong bond length distribution becomes 
consistent with a random graph dropped on top of a 
toroidal geometry. These bare lengths I should scale 
with the inverse square root of the remaining density, as 
they do, for example, at the TFIM fixed point. 

We take the 'long bond disease' exhibited by the bosons and 
Fibonacci anyons to indicate that the marginal fixed point be- 
havior they exhibit is actually caused by the planar truncation. 
Indeed, models which flow to weaker disorder generate many 
strong bonds at each step of the RG, of which many will be 
cut by the planar truncation rule, falsely preventing the flow 
back to weaker disorder. 

Finally, we note that the Ising anyon models have their own 
pathology within the planar approximation: between the op- 
posite sublattice rule and planarity, these models generate so 
few bonds that they fall apart into disjoint chunks at a finite 
time in the SDRG flow. This is also a physically spurious re- 
sult. Indeed, even the numerically observed narrowing of the 
bipartite Ising anyon coupling distribution without the planar 
approximation is believed to be a transient in what should be 
a marginal flow to infinite randomness.^ 



IV. MODELS 

For all of the models we consider, we use the following no- 
tational conventions. Latin coupling constants (J, h) refer to 
energetic coupling constants in the model and Q denotes the 
current maximum energy coupling. Greek letters (/?, Q refer 
to log-couplings relative to Q. Qq is the initial (unrenormal- 
ized) strongest coupling and thus T = log ft /ft measures the 
log energy scale of the flow. The number of sites remaining in 
the system is given by N. T,fl and the various coupling dis- 
tributions generally depend on the time in the renormalization 
flow t; we will generally omit writing this dependence except 
where necessary. 



A. Transverse Field Ising Model 

The strong disorder regime of the transverse field Ising 
model (TFIM) has been studied in great detail both an- 
alytically in one dimension 28 and numerically in higher 



a From Ref . 26 
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Hardcore bosons 


Planar 
Nonplanar 
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Steep Weakening 


TFIM 


Planar 
Nonplanar 
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0.43 

0.42 ± O.ofl 


Ising anyons 
(Nonbipartite) 


Planar 
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Nonexistent 
Weakening 


Ising anyons 
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Planar 
Nonplanar 
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(I) ~ iv- 1/2 


Nonexistent 
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Fibonacci anyons 


Planar 


Long bond disease 


"marginal" 





TABLE I. Behavior of various models under numerical SDRG in two dimensions. 



dimensions ! 26 I 29 E ] j n particular, Motrunich et clIf 6 demon- 
strated that in two dimensions, the TFIM flows to an infinite 
randomness critical point separating random quantum ferro- 
magnetic and paramagnetic phases - as is known analytically 
in one dimension. The TFIM is thus the fiduciary starting 
point for our treatment of 2-D models: it provides a check 
for both our numerical implementation and the validity of the 
planar approximation. If we are to have any faith in the pla- 
nar approximation in the more complicated topological mod- 
els below, it must reproduce the universal properties of the 
TFIM. 



1. Model 

The transverse field Ising model consists of a lattice of spin- 
1/2 degrees of freedom G{ governed by the Hamiltonian 



flipping of Gi : 



H 



■£■* 

(ij) 



ij^iVj 



^hiof-H^ma? 



(3) 



Here the J^ are Ising couplings, hi are local transverse fields 
and H is an externally applied uniform field in the z direction 
that couples to the moments /^ of the spins. In the strong 
disorder regime, the J^ and hi are broadly distributed random 
variables whose signs are irrelevant to the thermodynamics 
- any frustrated loop will have an extremely weak bond that 
may be dropped to unfrustrate the system (although clearly 
this effects the growth of net moments). In this regime it is 
natural to introduce the logarithmic variables: 

pi = log n/hi 
dj = log n/Jij 

where Q is the maximum coupling in the system at a given 
time in the flow. 

This system has two kinds of strong-disorder RG rules de- 
pending on whether the strongest coupling remaining in the 
system at scale T = \ogfto/ft is a transverse field hi or a 
bond Jij : 

• Field decimation: the field hi pins the spin Oi in the x 
direction so the site i may be removed from the lattice. 
For each pair j,k of o^'s former neighbors, perturba- 
tion theory generates an effective interaction by virtual 



J ik — Jjk 



J ji^ik 



jk " Jn ' 2hi 
or, in terms of the logarithmic variables: 

Cjk = min {(jk, Qi + Qk - A} 



(4) 



(5) 



• Bond decimation: the spins c^, dj coupled by the strong 
bond J^ bind together as a single larger Ising moment: 

Vi = Vi + H ( 6 ) 

This moment feels an effective transverse field 



h 



f r^i rij 



ZJi 



(7) 



and has couplings to all of the neighbors of either of the 
original spins 



J ik — ik * jk 

In terms of logarithmic variables: 

Pi = Pi + ft-Cij 

dk = min {dkXjk} 



(8) 



(9) 
(10) 



In the planar approximation, only the first rule needs to be 
modified; bond decimation simply collapses two sites joined 
by an edge and does not generate non-planar bonds. When a 
field pins a site with n neighbors, however, all of the (™) pos- 
sible bonds are in principle generated. In order to maintain 
local planarity, we therefore add these bonds in from strongest 
to weakest only when they would not cross a previously added 
bond. We note that Motrunich et al also threw out bonds that 
were generated by the RG rules by ignoring all bonds weaker 
than a certain threshold. This does not modify the results of 
the RG above the scale of the threshold but prevents an explo- 
sion in memory requirements. 




Flow of /3/C intercept and mean, TFIM 



FIG. 3. The facet flipping move used to 'mangle' initial lattices. By 
randomly introducing a finite density of these moves, we maintain 
planarity while broadening the degree distribution toward that ob- 
tained at a fixed point. 



2. Results 

In the numerical SDRG approach, one creates a bare two 
dimensional lattice G of some size Lx L with randomly sam- 
pled couplings j3 and ( from some Pq (/?) , Pq (£) and then runs 
the decimation procedure directly while monitoring the flow 
of geometry and couplings. Any infinite randomness fixed 
points will have much more complicated joint distributions 
Poo (G, /3, C) governing their couplings and geometry and thus 
the numerically observed flow will necessarily exhibit a tran- 
sient behavior as it approaches the scaling regime. Similarly, 
at the tail of the process as N becomes very small, we expect 
to see finite size effects modifying the flow. Thus, in order 
to identify thermodynamic scaling behavior, we would like to 
see scaling in the coupling distributions Pr(P), ^r(C) f° r as 
many orders of magnitude in N as possible between these two 
regimes, independently of the initial size L. 

All of our flows begin with periodic triangular lattices of 
size IxL, which we then 'mangle' randomly by applying 
facet flips to randomly chosen edges as in Figure [3] By ap- 
plying a finite density of these flips, we broaden the initial de- 
gree distribution of G (which for a triangular lattice is simply 
P{d) = S(d—6)) and bring it closer to the scaling distribution. 
Empirically, this suppresses the initial transients observed in 
the flow. 

The primary knobs for exploring the phase diagram of the 
TFIM are the initial distributions Ro((3) and Po(()- As in 
previous work, we find that the field distribution Rt(P) re- 
mains approximately exponential throughout all of our flows 
and thus we always take the initial condition 



Ro(P) = e 



-P 



(11) 



We also find that the bond distribution Pt(() tends to develop 
an upward initial slope as it flows towards the infinite random- 
ness critical point but this distribution gets cut off by a roughly 
exponential tail by the planar approximation. Thus, we take 



Po(C)^(a + 6C)e- cC 



(12) 



such that we can begin with an initial intercept a and slope 
m = b — ac. We locate the critical point by monitoring the 
intercept ratio R(0)/P(0) and the mean ratio (/3)/(C) anc * try 
to find fixed points in their flow as a function of initial condi- 
tions. See Figure [4] 

The primary features of the disordered TFIM in the planar 
approximation are consistent with those found by Motrunich 
etal.W 




FIG. 4. Rough cut through phase diagram near IRCP by varying 
intercept and slope (a, m) of linear exponential initial conditions for 
P(C) — ( a + K) e ~ c( * • The intercept ratio Ro/Po flows to large 
values in the paramagnet and small values in the ferromagnet and 
remains finite at the IRCR The mean ratio (f3) / (£) flows the opposite 
way. We identify our candidate IRCP at (a, m) ~ (0.10, 0.11) by 
looking for the separatrix in these flow lines. 
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FIG. 5. Flow of log field distribution R(/3) at putative IRCP. Ini- 
tial lattice 500 x 500, average over 400 runs, (inset) Flow of 
distribution width as measured by least squares fit to exponential 

R({3) = R oe - Rof3 and by flow of mean ({3}. 



1 . The model flows to a random ferromagnet, paramagnet 
or infinite randomness critical point depending on the 
relative strength of the initial bond and field distribu- 
tions. The distributions R(/3) and P(Q roughly scale at 
the IRCP. See Figs. [5] and § 

2. The infinite randomness critical point has a critical ex- 
ponent ip defined by w ~ N~^/ d where w is some 
measure of the width of the coupling distributions. 
Asymptotically all reasonable width measures should 
give the same exponent but they may be estimated nu- 
merically in various fashions. Using the same technique 
as Motrunich et. al, (least squares fitting to an expo- 
nential distribution for P(f3)), we find ip « 0.43, which 
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FIG. 6. Scaled flow of P(C/(C)) at putative IRCR Same dataset as 
Fig. [5] (inset) Flow of width as measured by mean {(). 
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FIG. 7. Flow of mean magnetic moment (/x) at putative IRCR Pro- 
vides estimate of fractal dimension df ~ 1.03 from (/x) ~ N~ d ^ d . 



agrees with Motrunich's result ip = 0.42 ± 0.06 and 
also with Monte Carlo and more recen t numerical re- 
finements of the SDRG calculation 2930 . Alternatively, 
one may estimate ijj by the flow of (() (or (/?)), which 
measures the width of the bond (field) distribution. This 
provides a higher estimate of i/j « 0.50 (0.49), still con- 
sistent with the earlier estimates to within error bars. 



3. The fractal dimension of the spin moments /x ~ 
jy-df/d ma y ^ e estimated from the flow of (/x) with 
A 7 '. We find df ~ 1.03 in agreement with other work. 
See Fig. [T] 



The planar SDRG accurately reproduces the phase diagram 
of the random TFIM as well as the critical exponents of the 
IRCP as determined by previous numerical studies. This in- 
dicates that the truncation of nonplanar bonds is irrelevant to 
the flow of a system with a strong infinite randomness fixed 
point. 



FIG. 8. Mean integrated out bond length during IRCP flow. This 
is a windowed time average of the geometric length of the bonds 
integrated out in a certain window of the flow. The horizontal line 
indicates the saturation line expected from a random graph on a torus 
of size L = 500; the dashed power law fit gives (I) ~ A -0-42 . 



B. Hardcore Bosons (XX Model) 

Disordered quantum XYZ chains have a rich phase diagram 
studied in considerable detail by Fisher. 28 The existence of in- 
finite randomness fixed points is a hallmark of these one di- 
mensional models. In higher dimensions, Motrunich et af^* 
did not find a strong numerical flow of the XX model (aka. 
hard-core bosons) to infinite randomness and did not pursue 
the model in much detail. We confirm this flow toward weaker 
disorder with a straightforward implementation of the boson 
SDRG. However, in the planar approximation, we find a spu- 
rious but numerically stable marginal fixed point at infinite 
randomness. This fixed point is characterized by saturation of 
the lengths of the strong bond and the breakdown of the bare 
geometry of the system. In this sense, the 'fixed point' behav- 
ior is not thermodynamic but rather that of a finite size system 
that has saturated. 



1. Model 

The random quantum XX model consists of a lattice of spin 
1/2 degrees of freedom c^ governed by the Hamiltonian 



(ij) 



(13) 



where J^ are random independently distributed couplings. 
Here we consider only positive couplings J^ and thus the 
ground state of any single term in the Hamiltonian is a sin- 
glet. 

There is only one strong-disorder RG rule for the XX 
model: 

• Bond decimation: the spins (Ji,(Jj form a singlet and 
drop out of the effective description of the system. For 
each pair fc, / of their former neighbors, perturbation 



theory generates an effective interaction by virtual ex- 
citations of this singlet: 



J'kl — Jkl 



JlJk 

Jij 



(14) 



where J\ ( J&) is the largest coupling between I (k) and 
either of the pair z, j. In log variables, 



Cm = min {Ckh (k + &- dj} • 



(15) 



The flow to an infinite randomness critical point of the 
TFIM arises when field and bond decimations exactly bal- 
ance throughout the RG flow - with only one RG rule, the XX 
model cannot exhibit such a behavior. Rather, if we begin the 
scale-free flow with various initial conditions, the only possi- 
ble behaviors are that the bond distribution gets wider with- 
out bound (a standard infinite randomness fixed point), nar- 
rower without bound (a finite randomness or pure fixed point) 
or finds a fixed distribution with width of order the starting 
width (a marginal infinite randomness fixed point). 



2. Results 

We initialize the planar SDRG with a periodic triangu- 
lar lattice of size L x L which we 'mangle' randomly (see 
Sec. |IV A 2| ) to reduce observed flow transients. The ini- 
tial couplings ( are chosen from an exponential distribution 
Po(C) = e~^ - several other initial distributions produced 
qualitatively similar results. The hardcore boson model flows 
in three stages (see Fig.|9]top): first, the system flows strongly 
to weaker disorder over about a decade in system size; second, 
the coupling distribution Pt(() stabilizes and the system ex- 
hibits an apparent marginal infinite randomness plateau; and, 
third, the flow begins drifting due to finite size effects at the 
smallest sizes. We have checked that each of these stages of 
the flow is consistent with the interpretation in terms of initial 
transient and finite size tail by simply varying the initial sys- 
tem size and comparing the onset and duration of each stage. 

The interpretation of the plateau during the second stage as 
evidence for a marginal infinite randomness fixed point must 
be physically spurious, as the system in the absence of the pla- 
nar truncation exhibits a consistent strong flow to weaker dis- 
order (data not shown). To understand the plateau better, we 
consider the geometry of the bare bonds during the flow and 
find that the mean bond length (/) (see Fig.[9]bottom) saturates 
the system size at the end of the first stage of the flow. After 
this point, the bond length distribution (not shown) is consis- 
tent with that of a random graph dropped onto a toroidal ge- 
ometry. In a physically correct fixed point, these bare lengths 
should scale with the inverse root of the density of remaining 
nodes, as they do at the IRCP of the TFIM (Fig. [8). This 'long 
bond disease' provides a diagnosis of the failure of the RG 
flow in the Fibonacci case as well. 
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FIG. 9. (top) Flow of mean coupling ( — InQ/J and (bottom) 
mean bond length for hardcore bosons in the planar approximation. 
Average of 20 planar SDRG runs for each initial lattice size (250 x 
250, 750 x 750), initial conditions as in text. The dashed horizontal 
lines indicate the expected distance between two randomly chosen 
points on a torus of size L. 



C. Ising Anyons 



Ising anyons are described by an SU(2) Chern- Simons the- 
ory at level k = 2. They may be equivalently be viewed as 
neutral Majorana fermion zero modes and their pairwise in- 
teractions may be recast as a free fermion hopping problem 
for real fermions. The undisordered Majorana chain describes 
a quantum critical point of a spinless superconductor 32 while 
the disordered model has the same strong-disorder rules as 
for hard-core boson hopping and therefore exhibits the same 
infinite randomness behavior as for the XX chain. 28 The two- 
dimensional bipartite case is equivalent to the bipartite imag- 
inary random hopping problem which is known to have a 
marginally stable infinite randomness phase 27 by various an- 
alytic mappings - this marginal flow is, however, hard to ob- 
serve in numerical SDRG investigations. We study this model 
on bipartite and non-bipartite lattices with and without the pla- 
nar approximation. 



1. Model 



The Ising any on model consists of a lattice of pinned any- 
onic degrees of freedom whose Hilbert space can be built up 
from the fusion rule 



CT(g)Cr = 10^ 



(16) 



where a represents a Ising any on and 1 , ip both represent vac- 
uum states (for our immediate purposes). The Hamiltonian 
that governs the model is then a sum of pairwise interactions 
that project a given pair onto one of its two fusion channels. 
Writing U.f- for the projector onto the i/j channel of fusion of 
the z, j pair of underlying particles, the Hamiltonian is 



tf = E J * 



IF 



(17) 



Since the 1 and tp fusion channels are both singlets which act 
trivially when fused with additional a particles, the sign of the 
interaction J^ is unimportant to the strong-disorder flow. 

There is only one strong-disorder RG rule for the Ising 
model: 

• Bond decimation: the anyons c^, (jj form a singlet and 
drop out of the effective description of the system. For 
each neighbor k of i and / of j , perturbation theory gen- 
erates an effective interaction by virtual excitations of 
this singlet: 



V 

J ki 



Ji 



ki 



JlJk 



(18) 



The rule is identical to that of the hard-core boson model ex- 
cept that bonds are only generated between next neighbors on 
opposite (local) sublattices. In particular, if the system begins 
with a bipartite lattice, the lattice remains bipartite through- 
out the flow and thus there are two distinct cases for the Ising 
any on RG: bipartite and non-bipartite. 



2. Results 

As expected, in the absence of the planar approximation, 
we find that both bipartite and non-bipartite Ising anyon mod- 
els flow slowly to weaker disorder, independent of initial con- 
ditions. No universal behavior can be extracted, but we in- 
dicate typical flows of the width as measured by the me- 
dian Ci such that 1 bond/site has strength greater than d 
in Fig. 
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These flows begin with either triangular (non- 
bipartite) or square (bipartite) lattices with couplings ( drawn 
from a linear-exponential distribution Po(C) — ( a + K) e ~ c ^ 
with intercept 1, slope at intercept 2. The shallow weakening 
of the bipartite Ising anyon distribution is consistent with the 
result that the actual thermodynamic flow is to a marginally 
stable IRFP, 27 which is too weakly attractive to be found at 
the sizes we consider. The non-bipartite flow is nearly as shal- 
low, which is suggestive that a similar behavior holds for the 
case of non-bipartite Ising anyons. This is not quite true, how- 
ever, as evinced by mappings of the non-bipartite pinned Ising 
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FIG. 10. (top) Representative flow of width as measured by median 
£i (1 bond/site has strength greater than £i) for bosons, non-bipartite 
Ising anyons, and bipartite Ising anyons without the planar approxi- 
mation. Averaged over 50 samples for each size (squares, L = 250, 
circles, L = 750) and model, (bottom) Mean integrated out bond 
length for same flows, showing saturation in the bosons. 



anyon problem into the disordered free fermion problem in 
symmetry class D, where a disordered gapless phase exists,^ 
but which does not exhibit infinite randomness scaling. 

Within the planar approximation, both bipartite and non- 
bipartite Ising anyon models 'fall apart' at a finite time (deci- 
mation fraction) in the RG flow. This arises due to the extreme 
truncation imposed by the opposite sublattice rule and the pla- 
nar approximation - a finite local sequence of decimations 
can leave behind disconnected sites. Thus, at a finite deci- 
mation fraction (of roughly 1.7%, measured for initial sizes 
L = 300, 500, averaged over 40 runs each), such motifs re- 
duce the system to a collection of disconnected sites and the 
RG stops prematurely. 



D. Fibonacci Anyons 

The Fibonacci anyons are described by an SU(2) Chern- 
Simons theory at level k = 3. They have attracted 
some interest in the context of topological quantum com- 
puting proposals 9 since they are the simplest non-Abelian 
anyons whose braiding rules are universal for quantum 
computation. 33 They are also candidate quasiparticles for the 
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FIG. 11. Strong-disorder rules for Fibonacci anyon models. In each configuration on the left, the thick bond Ji indicates the strong bond to be 
integrated out. For the chain, star and cross configurations, a strong (anti-)ferromagnetic bond leads to the configuration in the lower (upper) 
row with renormalized coupling J' F {J'a)- Positive couplings are antiferromagnetic (favor the vacuum channel). In the infinite disorder limit, 
the factor of the golden mean <j) — 1+ 2 ma y ^ e dropped, but the signs remain relevant. 



v = 12/5 fractional quantum Hall phase. 8 34 The Fibonacci 
chain, aka. the golden chain, was the first anyon lattice model 
to be studied in detail, 1214 both in its pure and disordered 
form, and exhibits an infinite randomness phase. In two di- 
mensions, this is the simplest model which suffers from the 



second difficulty described in Sec. [TlT| fusion of a pair of 
anyons may energetically influence the fusion of a disjoint 
pair of anyons elsewhere in the system (see the 'cross' config- 



uration in the rules, Fig. 11). Thus, this is the model which ex- 
hibits the full topological non-locality that motivated the use 
of the planar approximation in our treatment. 



1. Model 

The Fibonacci anyon model consists of a lattice of pinned 
anyonic degrees of freedom whose Hilbert space can be built 
up from the fusion rule 



TXT 



(19) 



where r represents a Fibonacci anyon and 1 represents the 
singlet or vacuum state. The Hamiltonian that governs the 



model is then a sum of pairwise interactions that project a 
given pair onto one of its two fusion channels. Writing H\- 
for the projector onto the r channel of fusion of the z, j pair 
of underlying particles, the Hamiltonian is 



h = j2j 



n r 



(20) 



Thus a positive coupling J^ corresponds to an "antiferromag- 
netic" term whose ground state is the "singlet" 1 fusion chan- 
nel. Restricted to a planar lattice of interactions, this is a suf- 
ficient description of the Hamiltonian. In general, however, 
the pair- wise interaction of two distant anyons should also be 
labeled by a path connecting the two anyons which indicates 
on which side of the other anyons in the system the interaction 
is mediated. 

The appropriate strong-disorder RG rules can be derived 
from an application of perturbation theory to various clusters 
of r particles. Unlike in the TFIM and XX models, the sign of 
the interaction is critically important to the physical behavior: 
a strongly coupled pair of rs can either fuse to a single new r 
or to the vacuum 1 state. In both cases new neighbor interac- 
tions are generated but with signs that depend on the topology 
of the local interaction network. Moreover, crossed bonds that 
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involve otherwise non-interacting r's interact under renormal- 
ization due to the non-locality of the Hilbert space; in partic- 
ular, decimation of a crossed bond flips the sign of the re- 
maining interaction and/or generates multiple path dependent 
interactions. The rules are summarized in Fig.fTT) 



2. Results 

As in the other models, we initialize the Fibonacci any on 
planar SDRG with a periodic triangular lattices of size L x L, 
which we 'mangle' with a finite density of facet flips (see 
Sec. |IVA"2| ). We choose independent random couplings ( for 
each bond from an exponential or linear-exponential distribu- 
tion Po(C) oc (a + b()e~ c ^. We tune these initial conditions 
in order to minimize the transient in the flow behavior and see 
as many decades as possible of scaling behavior. 

Under the planar SDRG, the Fibonacci anyons exhibit a 
spurious flow to a stable fixed point in which the coupling 
distribution is essentially fixed (full distribution not shown). 



See Figure 12 The plateau in the coupling width after an ini- 
tial transient is a robust feature of the planar flow, independent 
of the initial size and therefore suggestive of thermodynamic 
behavior -just as the upturn at small sizes shifts with the ini- 
tial size, which agrees with the interpretation of this upturn as 
a finite size effect. Unfortunately, the underlying geometry of 
this fixed point reflects the 'long bond disease' also exhibited 
by the boson flow and inconsistent with a physical fixed point 
(see Sec. |IVB2 ). As can be seen in Fig. [12J the mean bond 
length (I) of the strongest bond in the system saturates to the 
bare system size at essentially the same point in the flow that 
the coupling width plateaus. 



V. IMPLEMENTATION 
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FIG. 12. (top) The flow of the coupling width (£) downward to a 
plateau is falsely suggestive of a marginal infinite randomness fixed 
point, (bottom) The saturation of the mean decimated bond length 
(I) reflects the pathology of this fixed point, which has lost all sense 
of the two-dimensionality of the initial system. Average of 50 pla- 
nar SDRG runs for each initial lattice size (500 x 500, 750 x 750), 
initial conditions as in text. The dashed horizontal lines indicate the 
expected distance between two randomly chosen points on a torus of 
size L. 



The primary obstacle to implementing a planar strong- 
disorder renormalization group is the maintenance of the pla- 
narity condition on the underlying interaction graph. The 
usual representations of the interaction graph by either adja- 
cency matrix or edge lists does not contain information about 
the embedding of the graph into a two dimensional manifold, 
which quickly becomes impossible to discern. Similarly, us- 
ing a geometric embedding into a planar or toroidal surface 
is problematic because the bare geometry of the interaction 
graph may become very complicated during the renormaliza- 
tion flow. Rather, in order to maintain the local planarity of 
the interaction graph, we view it as the edge mesh on a topo- 
logical^ oriented surface represented using a halfedge data 
structure, familiar in the computational geometry literature. In 
this structure, interactions between vertices are represented by 
pairs of half edges which define the orientation of the adjacent 
facets on the surface (see Fig. 1 of Ref. 35). As the inter- 
actions and sites are decimated, we update this data structure 
in a topology preserving manner and thus guarantee that we 
always maintain local planarity correctly. 

We have implemented our own halfedge data structure li- 
brary in scientific Python 36 and run all planar renormaliza- 



tion groups using this code. The models studied without the 
planar approximation were also implemented in Python using 
the significantly simpler representation provided by adjacency 
lists. 



VI. CONCLUSIONS 

In summary, we have introduced the planar approxima- 
tion into the strong-disorder RG treatment of two dimensional 
pinned anyons in order to handle the crossing interactions gen- 
erated by the Fibonacci any on RG rules. By comparing the 
SDRG flow with and without the planar approximation on a 
variety of models, we find that the approximation accurately 
reproduces physics near systems flowing strongly to infinite 
randomness, such as at the random TFIM critical point, but 
breaks down in one of several characteristic ways in systems 
flowing back to weaker disorder. The similar pathologies of 
the hard-core bosons and the Fibonacci anyons strongly sug- 
gest that the Fibonacci model has no infinite randomness fixed 
points. 
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If a method can be developed to treat the topologically 
non-local interactions of the Fibonacci model without a pla- 
nar truncation, it is conceivable that it would reverse our con- 
jectured flow to weaker disorder. We believe this is unlikely 
and that the similarity to that of the 2-D hardcore bosons will 
continue to hold: the Fibonacci anyons do not exhibit any in- 
finite randomness physics. However, this does not rule out 
the possibility of other, nontrivial disorder-induced gapless 
phases in the Fibonacci phase diagram, as have been recently 
established in the form of a thermal metal phase for disordered 
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pinned Ising anyons. 
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